Report for Model 16

— for both jointly and separately

2024-11-08

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 4000000 iteration;
  • 3000000 burn-in;
  • 1000 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum      0.999       1.00
## LDevsum      1.001       1.01
## dh0          1.003       1.01
## dh1          1.002       1.01
## dh2          1.006       1.03
## dh3          1.003       1.00
## dl0          1.001       1.01
## dl1          1.003       1.02
## dl2          1.005       1.03
## dl3          1.003       1.02
## 
## Multivariate psrf
## 
## 1.01

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1208.16991 22190.0087
DIC3 1150.94079 22211.4544
PWAIC 41.54979 245.9433
WAIC 1176.99510 22233.4022

2. Separate Modeling of High-Level

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 4000000 iteration;
  • 3000000 burn-in;
  • 1000 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.00       1.00
## dh0           1.02       1.06
## dh1           1.04       1.15
## dh2           1.02       1.09
## dh3           1.00       1.00
## 
## Multivariate psrf
## 
## 1.04

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

H4
DIC 1358.50636
DIC3 1235.81461
PWAIC 86.58346
WAIC 1302.12921

3. Separate Modeling for Low-level

3.1 Trace plots for convergence check

The current MCMC setting is:

  • 4000000 iteration;
  • 3000000 burn-in;
  • 1000 thinning.

3.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## LDevsum      1.001      1.001
## dl0          1.004      1.022
## dl1          1.004      1.020
## dl2          0.999      0.999
## dl3          1.001      1.008
## 
## Multivariate psrf
## 
## 1.01

3.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

3.4 WAIC results

L7a
DIC 22335.6858
DIC3 22316.9560
PWAIC 315.3023
WAIC 22353.2432